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r^ , We extend, to include the effects of finite temperature, our earlier study of the interband dynamics 

On ' of electrons with Markoffian dephasing under the influence of uniform static electric flelds. We use 

^\ I a simple two-band tight-binding model and study the electric current response as a function of 

T-H . field strength and the model parameters. In addition to the Esaki-Tsu peak, near where the Bloch 

.^.j ' frequency equals the damping rate, we find current peaks near the Zener resonances, at equally 

O , spaced values of the inverse electric field. These become more prominenent and numerous with 

^^ ■ increasing bandwidth (in units of the temperature, with other parameters fixed). As expected, they 

^_. ' broaden with increasing damping (dephasing). 

'""' ! PACS numbers; 73.20.Dx, 71.70.Ej, 73.40. Gk 
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Q>^ ■ I. INTRODUCTION 

(n: 

^— ^ [ In a previous papem (henceforth referred to as I) we began a study of the effects on interband transitions of the 

1^^ • scattering of electrons from static imperfections in a semiconductor superlattice. A uniform static electric field was 

Qs ' applied. It has long been recognized that scattering destroys the coherence necessary to sustain the Bloch oscillations 

"T-*, , predicted in such a field, and, in practice this delayed the observation of these oscillations until the development of 

C^ ' semiconduckMr superlatticesoH. But the scattering destroys other types of interesting coherent motion, as well. Dunlap 

C , and KenkreHfl, in particular, looked at the effects on dynamic loeaUpation of electrons in time periodic electric fields. 

I ' We were interested in I in a multiband effect, the Rabi oscillationsBLl of electron population between bands of a crystal 

■ O in a uniform static electric field of appropriate magnitude — near an|^voided crossing of the two interpenetrating 

^ Wannier-Stark Ladders (WSL) arising essentially from different bandsQ. These are the so-called Zener resonances. 

Within a simple two-band tight-binding model we demonstrated the destruction of localization (of occupied electron 

states) to a single miniband by the dephasing associated with the scattering. The decay rate for the approach to steady 

state band populations exhibits sharp peaks at values of the static electric field which give Zener resonances in the 

absence of scattering. The specific approximations in I, however, effectively limited the results to the case of infinite 

^ ■ temperature. In particular, the steady state was assumed to be equal population of the two bands. In this paper 

Cu ' we remove that restriction, to discuss these effects at finite temperature. The necessary modification was described 

in Ref. ( ]^, namely relaxation of band populations toward values set by the Boltzmann factors describing thermal 

equilibrium. This will allow us, in particular, to look at the electric current, a quantity of obvious experimental 

interest which vanishes in the infinite temperature limit. 

II. MODEL 

We consider the same model that we treated in I, a standard simple tight-binding modelcl of a two-band system in 
a static electric field E. The Hamiltonian can be written as 

i/ = ^ (Aa + nujB)al^an + (A^ + nujB)b\jjn 

-iWa/A){ai+^an + h.c.) + (VKfc/4)(6,\+i6„ + h.c.) 
+eER{albn + blun) 



o 



>< 



(1) 



Here the subscripts label the lattice sites and the lower and upper minibands are designated by symbols a and 6, 
respectively. We have introduced the notation ujb = eEd, where d is the lattice constant, for the Bloch frequency, 
which will appear often below. The first two terms describe the site energies of the Wannier states in the presence of 
the electric field, and Wa,b are the widths of the isolated {E = 0) minibands induced by nearest neighbor hopping: 
£a,bik) — Aa^b T (W^a,&/2) cos/c, whcrc the dinicnsionless wave vector k is in units of the inverse lattice constant d. 
The last term is the on-site electric dipole coupling between minibands; eR is the corresponding dipole moment. This 
Hamiltonian does neglect Coulomb interactions and electric dipole elements between Wannier states on different sites, 
but it contains the essential physics for the problemQ^Ej. Note that the hopping parameters Wa,b are written here 
with opposite signs, so that with both parameters positive the band structure at _E = is of the standard nearly free 
electron character, with direct band gaps at the zone boundary. But the calculation to follow is valid for arbitrary 
signs of the parameters. 

It is easily shownEl'ty that the exact spectrum of H is two interpenetrating Wannier-Stark Ladders. But what do 
the corresponding states represent in terms of the occupation of the original bands as a function of time, and what is 
the influence of scattering? For vanishing dipole matrix element between bands, R = 0, there is no interband mixing. 
Each of the two bands gives rise to a single WSL. Clearly, when the electric field amplitude is such that the ladders 
become degenerate, even small values of R lead to strong interband mixing. The crossing of the ladders is "avoided" 
by any finite R, and the behavior near those avoided crossings (the Zener resonances) is of particular importance and 
interest. In general there are peaks in the current response at those values of the electric field, but the peaks are 
broadened by increasing temperature, as well as by decreasing bandwidth relative to band separation. 

We start by defining the density matrix in the representation of the two bands. 

Pit) = E P^nCa. (2) 

ijmn 

where i,j — 1 or 2 are band indices: ^^ {^j^^) and ^^ (^^) designate aj„ {am) and 53„ (bm), respectively. Since we 
are interested in the dynamics of occupation of various band states, it is convenient to work in a wave vector basis, 
by Fourier transforming the density matrix. In general, since pmn is not translationally invariant (a function only of 
m — n), we have a full set p'l — X^mn Pmnit) ex.p[—ikm + iqn] of Fourier components. But we will be interested in 

the wave vector diagonal band occupation numbers p]^i^{t) = p'-' (fc, t). Then at finite temperature T we insist that the 
wave vector and band diagonal occupation number relax to the thermal equilibrium value, 

pi?,^p^(fc)^e-/'^'W/Ee-''^^('\ (3) 

where f3 — l/kgT , the sum in the partition function is over j = 1,2 and over all wave vectors q, and the band energies 
are those given above, 

ei,2ik)^Aa,bTiWa,b/2)cosk . (4) 

Within a constant relaxation rate approximation^, the density matrix p{k, t) satisfies the following stochastic Liou- 
ville equation (SLE) (we set ?i = 1 throughout this paper), 

i^^[H,pit)]-^^[p{t)~pT]. (5) 

Here each of H, p, and pT is labelled by (the same) wave vector k. The operator F describes the relaxation of the 
off-diagonal elements of p through dephasing: 

r[p - Pt] = E"'^- ^''(^'*) - ^^jPTik)] eik)eik). (6) 

The utility of this simplest form of the SLE has been discussed by Kenkre and collaborators (see Ref. and references 
therein). The parameters a^ measure the loss of phase coherence between sites, or the scattering lifetime of band 
states labeled by quasimomentum. 

As in I it is convenient to introduce the linear combinations of density matrix elements: 

p+{k,t)=p'\k,t)+p''{k,t), (7) 

p-{k,t)=p'\k,t)-p''{k,t), (8) 

p+-ik,t)=p'\k,t)+p'\k,t), (9) 

p-+{k,t)=i[p'\k,t)~p'\k,t)]. (10) 



For simplicity, we also take an = a22 = 0^12 = C(2i = a to reduce the number of parameters in the theory. Then the 
SLE (g) has the explicit components 

d d 

—p+{k,t)-uiB-g^p+{k,t) = a[p+{k,t)- p^{k)] , (11) 

—p_ik,t)~iJBgi^P-ik,t) = -2eERp_+ik,t)-a[p_ik,t)-p^{k)] , (12) 



dt 



p+_ (fc, t) - ^BT^P+- (fc, t) - (a - H/ cos fc) p_+ {k,t)~ ap+^ {k,t) , (13) 



^^p_+{k,t) - LjBgj:P-+{k,t) ^ ^ - (^A-Wcosky+^ik,t) + 2eERp_{k,t)-ap_+ik,t) . (14) 

Here we have used the simplified notation A = Aq — A;, and W = {Wa + Wb)/2. The equation ( [ll|) for p+ is decoupled 
from the others, and is readily integrated to give 

p+{k,t)^e-°'Up+{k + uJBt)+a f dt'e°'''p+[k + u;Bit-t')]\. (15) 







The equations for p^{L.t), p^ (fc,i), and p i-(^i^) can be reduced to the following ordinary differential equations 

in an accelerated basislld, k{t) — k — ojBt or, equivalently, in the transverse or vector gauge discussed in Ref. |l^, 

—X(k,t) = -~2eER Z{k,t) -~a[X{k,t) - p^(/c - w^t)], (16) 

— y (fc, t) = [A-W cos(fc - ujBt)]Z{k, t) - aY{k, t), (17) 

—Z{k,t) = -[A - W cos{k - ujBt)]Y{k,t) + 2eER X{k,t) - aZ{k,t) . (18) 

Here X{k,t) = p-{k ~ LJBt,t), Y{k,t) = p^ [k — uJBt,t), and Z{k,t) = p ^{k — w^i, i). The structure is exactly 

the same as Eqs. (17 - 19) in I, except for the fc-dependent relaxation of X here. As we did there, we can integrate 
these equations analytically as a perturbation series in the parameter p = 2eER, which characterizes the electric 
dipole coupling between bands. To lowest nontrivial (second) order in p we find 

p^{k,t)^e~"^lp;j.{k + LJBt)+a dt' e"*' p;^[k + uJB{t - t')] - p^ dt' dt" pj.{k + ujBt) 
I Jo Jo Jo '- 



+ a dr e"'^pj, [k + t^s (t - t)] 

"'0 



cos I' 



dr' [A - -VKcos {k + ujB{t - t'))] \ . (19) 



III. CURRENT 

We turn now to the calculation of the interesting physical quantity, the current j{t) along the superlattice direction. 
In the band i = a,b the instantaneous current is given by the sum over wave vectors of the relevant electron velocity 
Vi{k) = (Wid/2) sink times the number density in that band at that wave vector, nQPi{k,t), where no is the number 
of carriers per unit area in each cell of the superlattice. In terms of the convenient quantities p±{k,t) we then have 

j{t)= —[{Wa-Wb)p+{k,t) + 2Wp^{k,t)] nod sink. (20) 

Jo Stt 

Of particular interest is the steady state long time average of this, 



0-)= lim ^ f "dtj{t). (21) 

To^oo io Jo 

To second order in /i we have the probabiHty densities p±{k, t), in Eqs. (Qq) and (US), and we find 



nod \ 2 J \a'^ + U!^ 
'W\ f an 



4 ; V "2 + wl 



[Cl + {Wb/Wa)C2] 

; oo 

[C1+C2] ^ [i?7+i?+] J,2(VF/a;s), (22) 



with 



^ e^/^A/i(/?M/j2) 

'-e-'3A/o(/3Tyj2) + /o(/3Ty,/2) ' ^ ' 

^ h[(iWb/2) 

^'^ - e-P^hiPWa/l) + h{PWb/2) ' ^ ' 

where /„ is the modified Bessel function and J„ the ordinary Bessel function of order n. The interband effects are ah 
contained in the second term on the right hand side of (|22|), proportional to /i^. The Zener resonances, near A = uujb, 
with n an integer, exhibit themselves as peaks in the factors D^ . 

We will look at the nonlinear conductance, the current as a function of increasing electric field, which is conveniently 
parameterized in the dimensionless form ws/A. At sufficiently high fields, as the Wannier-Stark functions become 
increasingly localized in space, the current falls off inversely with cob- To the extent that the interband contribution can 
be neglected (small /i) the field dependence of the response is given by the factor ^^/(a^ +ujg), with an "Esaki-Tsu" 
peakB near ujb — ct. 

There are two sources of temperature dependence in (|22|). The first is the explicit appearance of /3 = 1/kBT in the 
factors Cl and C2. The other is implicit; the relaxation rate a is ordinarily also temperature dependent. Though that 
can be modeled for a specific relaxation mechanism, we simply take it to be constant below. At high temperatures, 
with /3A, (3Wi <C 1, we have Ci and C2 approximately linear in (3\ the current (j) also therefore falls off as l/T. 
This is the well known effect of approaching uniform thermal population of the states throughout each band as the 
temperature rises; at complete uniformity the current vanishes. 

All of these effects are seen in Fig. 1, where the interband matrix element /x has been set equal to zero, and 
a/A = 0.1. The second term on the right hand side of Eq. (p3) is then absent, and the single remaining peak arises 
from the first term in that equation. To avoid a multiplicity of parameters we have taken the two bandwidths to 
be equal, Wa = Wb = W = 0.8A. In the remaining figures we have set the interband coupling to be /it/A = 0.4. 
Again Wa — Wb — W, with increasing values W/A — 0.4 and 0.8 in Figs. 2 and 3, respectively. Each figure gives 
results for three values of the parameter /3A, namely 1, 10 and 100. The value of unity is representative of a typical 
superlattice, of lattice parameter of order 100 A and carrier effective mass m* /m « 0.1, at room temperature. From 
Eq. ( p2|) we see that the height in (j) of the nth peak (for values of n up to 5 or 6, where a^ ^ '-^b)^ ^t ujb — A/n 
is approximately proportional to [n/(l + .Oln^)]'^ J^_^{nW/ A) . The Bessel function overlap factors increase rapidly 
with increasing bandwidth W over the range we have chosen for W , and so we see more well-defined Zener resonance 
peaks with increasing bandwidth , going from Fig. 2 to Fig. 3. The major effect of increasing the damping rate a is 
to broaden the peaks in the current. Thus, for the particular choice of W/A = 0.4 and fi/A = 0.4, we can compare 
the results for a/ A = 0.05 in Fig. 4 and a/ A = 0.2 in Fig. 5 with the previous curves for a/ A = 0.1 in Fig. 2. 

IV. CONCLUSIONS 

Within a simple two-band tight binding model we have studied the electric current response of a semiconductor 
superlattice subject to a finite uniform electric field. Relaxation processes have been assumed Markoffian, and they 
have been described by a single parameter characterizing the rate at which a nonequilibrium density is restored to 
thermal equilibrium, using a stochastic Liouville equation. We have been particularly interested in the interband 



transitions, which lead to current peaks near the values of external field where Zener resonances occur. We have 
exhibited the variation in height and width of these peaks, as well as in the overall current magnitude with temperature 
and with the various parameters of the theory, including bandwidth, band separation, interband dipole coupling, and 
relaxation rate. 

The results we have obtained are perturbational in /^; we do require weak interband coupling. Moreover, we have 
limited the discussion to a single pair of bands, assuming the impact of all other bands to be negligible on the Zener 
tunneling between these two. This can be realized in practice by using a dimerized semiconductor superlattice with, 
for example, uniform wells but alternating thick and thin barriers between them. With weak alternation of barrier 
thickness one can adjust a pair of bands resulting from the doubling of the unit cell to be well isolated from all other 
bands, and the predictions of this paper can be studied in such a system by varying the external uniform electric 
field. Another way to observe the effects predicted here is to use ultracold atoms in accelerating optical potentials. 
Recently, Rabi oscillations were observed in such a system, where the interband coupling was generated by a small 
phase modulationt3. Very low temperatures can be realized in these systems, and it should be possible to observe the 
resonant structures in current that we have predicted in this paper. 

This work was supported in part by the National Natural Science Foundations of China under Grant No. 19725417, 
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FIG. 1. Time averaged current in the absence of interband transitions {fi = 0). The ratio of relaxation rate to band center 
separation is a/A — 0.1. The bandwidths are Wa = Wb = W = 0.8A. 

FIG. 2. Time averaged current, as in Fig. 1, but with fj, — W = 0.4A. 

FIG. 3. Time averaged current, as in Fig. 2, but with W — 0.8A. 

FIG. 4. Time averaged current, as in Fig. 2, but with a — 0.05A. 

FIG. 5. Time averaged current, as in Fig. 2, but with a = 0.2A. 
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